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The Ca+, Sr+, Ba+, and Yb+ ions immersed in an ultracold gas of the Cr atoms are proposed 
as experimentally feasible heteronuclear systems in which ion-atom interactions at ultralow tem¬ 
peratures can be controlled with magnetically tunable Feshbach resonances without charge transfer 
and radiative losses. Ab initio techniques are applied to investigate electronic-ground-state proper¬ 
ties of the (CaCr) + , (SrCr)'*', (BaCr) + , and (YbCr)'*' molecular ions. The potential energy curves, 
permanent electric dipole moments, and static electric dipole polarizabilities are computed. The 
spin restricted open-shell coupled cluster method restricted to single, double, and noniterative triple 
excitations, RCCSD(T), and the multireference configuration interaction method restricted to single 
and double excitations, MRCISD, are employed. The scalar relativistic effects are included within 
the small-core energy-consistent pseudopotentials. The leading long-range induction and dispersion 
interaction coefficients are also reported. Finally, magnetic Feshbach resonances between the Ca"*", 

Sr+, Ba+, and Yb+ ions interacting with the Cr atoms are analyzed. The present proposal opens 
the way towards robust quantum simulations and computations with ultracold ion-atom systems 
free of radiative charge transfer losses. 

PACS numbers: 34.20.-b,33.15.Kr,34.50.Cx,34.70.-|-e 


I. INTRODUCTION 

Hybrid systems of laser-cooled trapped ions and ultra¬ 
cold atoms combined in a single experimental setup are 
currently emerging as a new platform for fundamental 
research in quantum physics [1]. Potential applications 
range from controlled ion-atom collisions and chemical 
reactions [2-5] to quantum computations [6] and simu¬ 
lations of solid-state physics [7, 8]. Another possible re¬ 
search direction is the formation of cold molecular ions [9] 
for precision measurements [10], ultracold chemistry [11], 
or other novel applications [1]. 

Losses due to radiative inelastic processes, that is the 
radiative charge transfer and radiative association [4, 12], 
are the basis of controlled ion-atom chemistry [11] but 
can obstruct sympathetic cooling [3] , and interesting ap¬ 
plications in quantum computations [6] and quantum 
simulations of many-body physics [8]. Radiative losses 
can be reduced by choosing a system with a favorable 
electronic structure such as Yb+d-Li [9] and avoided if, 
and only if, an interacting ion-atom system is in its ab¬ 
solute electronic ground state and the charge transfer is 
impossible. This condition is met when the ionization 
potential of an atom is larger than the electron affinity 
of a positive ion. 

Up to date, all experimentally investigated cold het¬ 
eronuclear hybrid ion-atom systems have presented ra¬ 
diative losses [2-5, 13-21]. These losses are unavoidable 
for almost all alkaline-earth-metal and ytterbium ions im¬ 
mersed in a cold gas of alkali-metal atoms [9, 22-27] (ex¬ 
cept the Ba+ ion in a gas of the Li atoms). They can 
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be, however, avoided for some combinations of alkali- 
metal ions interacting with alkali-metal atoms [28] or 
alkaline-earth-metal ions interacting with alkaline-earth- 
metal atoms [18] and for almost all combinations of alkali- 
metal ions interacting with alkaline-earth-metal atoms 
(except the Li+ ion in a gas of the Ba atoms) [cf. Ta¬ 
ble I]. Homonuclear ion-atom systems do not suffer from 
radiative losses but show resonant spin-changing charge- 
exchange inelastic collisions [29-31]. Unfortunately, the 
preparation and manipulation of laser-cooled closed-shell 
alkali-metal ions and closed-shell alkaline-earth-metal 
atoms are experimentally more challenging than using 
open-shell ones. Furthermore, the existence of magneti¬ 
cally tunable Feshbach resonances predicted for mixtures 
of closed-shell and open-shell species is still waiting for an 
experimental confirmation [32, 33]. Therefore, the useful 
control of ultracold ion-atom interactions with magnet¬ 
ically tunable Feshbach resonances should be expected 
only when both ions and atoms are of the open-shell na¬ 
ture [34]. 

In the present work, we propose the Ca+, Sr+, Ba’*', 
and Yb+ ions immersed in an ultracold gas of the Cr 
atoms as systems in which ion-atom interactions at ul¬ 
tralow temperatures can be controlled with magnetically 
tunable Feshbach resonances without charge transfer and 
radiative losses. Among the present experimentally ac¬ 
cessible laser-cooled trapped ions and ultracold atoms, 
the other possible combinations are the Ba’*' ions and 
the Li atoms as well as the Sr+ or Ba’*' ions and the 
Dy or Er atoms. The interaction of the state alkaline- 
earth-metal or ytterbium ion with the state chromium 
atom gives rise to the two molecular electronic states 
of the A®S+ and a®S+ symmetries [35]. Here we re¬ 
port the ab initio properties of these electronic states of 
the (CaCr)+, (SrCr)+, (BaCr)+, and (YbCr)+ molecular 
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TABLE 1. The hybrid heteronuclear ion-atom systems without charge transfer radiative losses among the combinations of the 
present experimentally accessible laser-cooled trapped ions and ultracold atoms. 


Nature of ion/atom 
open-shell / open-shell 
open-shell / closed-shell 
closed-shell / open-shell 

closed-shell / closed-shell 


Systems without charge transfer radiative losses 

Ba+/Li, Ca+/Cr, Sr+/Cr, Ba+/Cr, Yb+/Cr, Sr+/Dy, Ba+/Dy, Sr+/Er, Ba+/Er 
Ca+/Mg, Sr+/Mg, Ba+/Mg, Yb+/Mg, Sr+/Ca, Ba+/Ca, Ba+/Sr, Ca+/Yb, Sr+/Yb, Ba+/Yb 
Na+/Li, K+/Li, Rb+/Li, Cs+/Li, K+/Na, Rb+/Na, Cs+/Na, Rb+/K, Cs+/K, Cs+/Rb 
alkali-metal ion/Cr, alkali-metal ion/Dy, alkali-metal ion/Er 
alkali-metal ion/alkaline-earth-metal atom (except Li’^/Ba) 


ions and prospects for magnetically tunable Feshbach res¬ 
onances opening the way towards more elaborate studies 
of the formation and application of the proposed hybrid 
ion-atom systems. 

The plan of our paper is as follows. Section II describes 
the theoretical methods used in the electronic structure 
and scattering calculations. Section III discusses the po¬ 
tential energy curves and electric properties of the stud¬ 
ied molecular ions and prospect for magnetically tunable 
Feshbach resonances. Section IV summarizes our paper. 


II. COMPUTATIONAL DETAILS 
A. Electronic structure calculations 

We adopt the computational scheme successfully ap¬ 
plied to the ground states of the (LiYb)+ molec¬ 
ular ion [9], the chromium-alkaline-earth-metal-atom 
molecules [36], and the europium-alkali-metal-atom 
molecules [37]. The potential-energy curves for the a®E+ 
electronic state in the Born-Oppenheimer approximation 
are obtained with the spin-restricted open-shell coupled 
cluster method restricted to single, double, and nonitera¬ 
tive triple excitations, starting from the restricted open- 
shell Hartree-Fock (ROHF) orbitals, RCCSD(T) [38]. 
The energy splittings between the potential energy curves 
of the and a®E+ electronic states are calcu¬ 

lated with the multireference configuration interaction 
method restricted to single and double excitations (MR- 
CISD) [39]. The potential energy curves of the A^E’*' 
electronic ground state are obtained by subtracting the 
energy splittings from the potential energy curves of the 
a®E+ electronic state calculated with the RCCSD(T) 
method. The interaction energies are obtained with the 
supermolecule method with the basis set superposition 
error corrected by using the counterpoise correction [40] 

^(XCr)+ = E(^xCr)+ — Ex+ - Eqt , (1) 

where E(^xCr)+ denotes the total energy of the molecular 
ion, and Ex+ and Ecr are the total energies of the ion 
and atom computed in the diatom basis set. 

The scalar relativistic effects are included by employing 
the small-core relativistic energy-consistent pseudopoten¬ 
tials (ECP) to replace the inner-shells electrons [41]. The 


use of the pseudopotentials allows one to use larger ba¬ 
sis sets to describe the valence electrons and models 
the inner-shells electrons density as accurately as the 
high quality atomic calculation used to fit the pseudopo¬ 
tentials. The pseudopotentials from the Stuttgart li¬ 
brary are employed in all calculations. The Cr atom 
is described by the ECPIOMDF pseudopotential [42] 
and the [14sl3pl0(i5/453/i] basis set with coefficients 
taken from the aug-cc-pVQZ-DK basis [43]. The Ca, Sr, 
Ba, and Yb atoms are described with the ECPIOMDF, 
ECP28MDF, ECP46MDF, and ECP28MDF pseudopo¬ 
tentials [44, 45] and the [12sl2p7dAf 2g], [14sllp6(i5/4g], 
[13sl2p6d5/45r], and [15sl4pl2dll/8g] basis sets, respec¬ 
tively, obtained by augmenting the basis sets suggested in 
Refs. [44, 45]. The basis sets are additionally augmented 
in all calculations by the set of [3s3p2dlflg] bond func¬ 
tions [46]. 

The interaction potential between the A+(^S') ion and 
the Ct(JS) atom, both in the electronic ground state, at 
large internuclear distances R reads [22] 

^ind /^ind ^disp 

= + ( 2 ) 

where the leading long-range induction coefficients are 
= ^acr and = ^/Scr, and the dispersion one is 

Cg*’*’’ = f /q°° Q;x+(* a;)acr(*a;)da;. acr is the static elec¬ 
tric dipole polarizability of the Cr atom, /3cr is the static 
electric quadrupole polarizability of the Cr atom, and 
ax+ is the dynamic polarizbility of the X~^ ion 

and the Cr atom at imaginary frequency. The dynamic 
polarizabilities at imaginary frequency of the atom and 
ions are obtained from the sum over state expression us¬ 
ing the transition moments from Refs. [36, 47, 48] and en¬ 
ergy levels from the NIST Atomic Spectra Database [49] . 
The static quadrupole polarizability of the Cr atom is 
calculated with the RCCSD(T) and finite field methods. 

The permanent electric dipole moments and static 
electric dipole polarizabilities are calculated with the fi¬ 
nite field method using the RCCSD(T) and MRCISD 
methods for the a®E+ and A®E+ electronic states, re¬ 
spectively. The z axis is chosen along the internuclear 
axis, oriented from the alkaline-earth-metal ion to the 
chromium atom, and the origin is set in the center of 
mass. 

All electronic structure calculations are performed 
with the Molpro package of ab initio programs [50]. 
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B. Scattering calculations 

The Hamiltonian describing the relative nuclear mo¬ 
tion of the X^(^S) ion and the Cr(^5') atom reads 

I 

^ = ^ ysmS,Ms){S,Ms\ + V^%R) + Hx+ + Her , (3) 

^ ^ SMs 


where /r is the reduced mass, R is the internuclear dis¬ 
tance, I is the rotational angular momentum operator, 
and Vs{R) is the potential energy curve for the state 
with total electronic spin S. is the spin-dipolc- 

spin-dipole interaction responsible for the dipolar relax¬ 
ation [51]. The atomic Hamiltonian, Hj {j = Cr), 
including Zeeman and hyperfine interactions, is given by 

Hj — Cj^j * Rz ( 4 ) 

where Sj and ij are the electron and nuclear spin op¬ 
erators, C.j is the hyperhne coupling constant, g^/j are 
the electron and nuclear g factors, and /i^ /n ^'I'e the 
Bohr and nuclear magnetons. The experimental val¬ 
ues of the nuclear g factors from the NIST Atomic 
Spectra Database [49] and the following hyperfine cou¬ 
pling constant are assumed: —806.4 MHz for ^^Ca+ [52], 
-1000.5MHz for ®^Sr+ [53], 4018.8MHz for [54], 

-3497.2 MHz for [55], and -83.60 MHz for 

®^Cr [56]. Negative values mean that the total angu¬ 
lar momentum is larger for the ground hyperfine state as 
compared to the excited one. The impact of the Lorentz 
force and Landau quantization effects [57] is neglected in 
the present study. 

The total scattering wave function is constructed in a 
fully uncoupled basis set assuming the projection of the 
total angular momentum to be conserved. The coupled 
channels equations are solved using the renormalized Nu- 
merov propagator [58] with step-size doubling. The wave 
function ratio is propagated to large interatomic separa¬ 
tions, transformed to the diagonal basis, and the K and S 
matrices are extracted by imposing the long-range scat¬ 
tering boundary conditions in terms of the Bessel func¬ 
tions. The scattering length is obtained from the S ma¬ 
trix for the entrance channel ao = (l — Sii)(l + Sii)/(ik), 
where k = y^2fiE/h? with the collision energy E [59-61]. 

III. NUMERICAL RESULTS AND DISCUSSION 
A. Potential energy curves 

The interaction of the ground-state alkaline-earth- 
metal or ytterbium ion with the ground-state chromium 
atom, which both are open-shell, results in the two elec¬ 
tronic states of the a®E+ and symmetries. The 
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FIG. 1. (Color online) Potential energy curves of the (a) 
and (b) electronic states of the (CaCr)"*", (SrCr)^, 

(BaCr)'*', and (YbCr)'*' molecular ions. 


computed potential energy curves for these electronic 
states of the (CaCr)+, (SrCr)+, (BaCr)+, and (YbCr)+ 
molecular ions are presented in Fig. 1(a) and Fig. 1(b), 
respectively, and show a smooth behavior with well- 
dehned minima. Their equilibrium distances Rg and well 
depths De are reported in Table H. The energy splittings 
between the a®E+ and states are shown in Fig. 2. 

The leading long-range induction and dispersion interac¬ 
tion coefficients are reported in Table HI. 

The ion-atom interaction is dominated by the induc¬ 
tion contribution, that results in relatively large bind¬ 
ing energy of the electronic ground state. The 

well depths of this state are between 7051 cm“^ for the 
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TABLE II. Spectroscopic characteristics: Equilibrium bond length R^, well depth De, harmonic frequency ojo, and number 
of bound vibrational states of the and a®E+ electronic states and rotational constant Bq, electric dipole moment 

do, average polarizability qq, and polarizability anisotropy Aao, for the rovibrational ground level of the and 

electronic states of the (CaCr)"*", (SrCr)"*", (BaCr)"*', and (YbCr)^ molecular ions. Masses of the most abundant isotopes are 
assumed. 


System 

Re (bohr) 

De (cm ^) 

cJo (cm ^) 

Ne 

Bo (cm ^) 

Irfol(D) 

ao (bohr^) 

Aao (bohr®) 





a®S+ 





(CaCr) + 

6.72 

3582 

87.7 

110 

0.059 

3.58 

158 

88 

(SrCr) + 

7.04 

3007 

66.5 

125 

0.041 

1.17 

176 

80 

(BaCr) + 

7.06 

3235 

67.2 

136 

0.032 

0.48 

187 

20 

(YbCr) + 

6.93 

2998 

58.9 

139 

X®E+ 

0.031 

1.18 

142 

81 

(CaCr)+ 

5.84 

10233 

147.5 

154 

0.078 

3.17 

126 

80 

(SrCr)+ 

6.13 

8534 

115.5 

171 

0.049 

1.22 

145 

104 

(BaCr) + 

6.45 

7051 

98.3 

172 

0.038 

0.45 

185 

157 

(YbCr) + 

5.91 

9879 

110.7 

202 

0.043 

2.19 

104 

60 


TABLE III. Induction and dispersion interaction coefficients 
describing the long-range part of the interaction potential be¬ 
tween the alkaline-earth-metal and Yb ions, and the Cr atom, 
all in the electronic ground state. 


System 

(a.u.) 

cr® (a.u.) 

Ct^ (a.u.) 

Ca+-bCr 

40.8 

290 

557 

Sr+-tCr 

40.8 

290 

654 

Ba+-tCr 

40.8 

290 

807 

Yb+-bCr 

40.8 

290 

484 



R (bohr) 


FIG. 2. (Color online) Energy splittings between the X®E+ 
and a®E^ electronic states of the (CaCr)"*', (SrCr)"*', (BaCr)"*', 
and (YbCr)'*' molecular ions. The points indicate the values 
for the equilibrium distance of the X^E"*"state. 


(BaCr)+ molecular ion and 10233 cm“^ for the (CaCr)+ 
molecular ion. The corresponding equilibrium distances 
are between 5.84 bohr and 6.45 bohr for the (CaCr)+ 
and (BaCr)’*' molecular ions, respectively. The bind¬ 


ing energy of the a®S+ high-spin electronic state is 
around three times smaller and its equilibrium distance is 
around 1 bohr longer as compared to the electronic 

ground state. The well depths of this state are between 
2998 cm“^ for the (YbCr)+ molecular ion and 3582 cm“^ 
for the (CaCr)'*' molecular ion. The corresponding equi¬ 
librium distances are between 6.72 bohr and 7.06 bohr 
for the (CaCr)+ and (BaCr)+ molecular ions, respec¬ 
tively. The characteristics of the presented potential en¬ 
ergy curves are similar with a small variation of the bind¬ 
ing energy and equilibrium distance among investigated 
ion-atom systems. The long-range parts of the interac¬ 
tion potentials are also almost identical [cf. Fig. 1 and 
Table III]. This behavior results from the fact that the 
electronic properties of the investigated ion-atom systems 
are determined by the induction interaction between dif¬ 
ferent ions (in the same electronic state) but the same 
chromium atom. 

The electronic ground state of the investigated 

molecular ions is bound around three times stronger 
as compared to the isoelectronic chromiuin-alkali-metal- 
atom molecules [62, 63]. The stabilization due to the 
induction interaction is even more pronounced in the 
case of the a®S+ high-spin electronic state which is 
bound around twenty times stronger in molecular ions 
as compared to their isoelectronic nentral counterparts. 
In fact, the potential energy curves for the a®S+ high- 
spin electronic state of the investigated molecular ions 
are as deep as for the X®S+ electronic ground state 
of the chromiuin-alkali-metal-atom molecules and the 
electronic ground state of the chromium-alkaline- 
earth-metal-atom molecules. Interestingly, the energy 
splittings between the and a®S+ electronic states 

(dominated by the exchange energy and responsible for 
standard magnetic Feshbach resonances) of the investi¬ 
gated molecular ions are of the same order of magni¬ 
tude as in the isoelectronic chromium-alkali-metal-atom 
molecules. 
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FIG. 3. (Color online) Permanent electric dipole moment of 
the (CaCr)”'", (SrCr)^, (BaCr)”'", and (YbCr)^ molecular ions 
in the (a) and (b) electronic states as a function 

of the internuclear distance. The z axis is oriented from X'^ 
to Cr and the origin is in the center of mass. The points 
indicate the values for the ground rovibrational level, and 
the dotted lines represent the permanent dipole moment that 
the molecular ion would have if the charge was completely 
localized on the ion. 


The accuracy of ab initio calculations depends on the 
proper treatment of relativistic effects, the reproduction 
of the correlation energy, and the convergence in the 
size of the basis function set. The employed CCSD(T) 
method is the golden standard of quantum chemistry 
and a good compromise between the computational cost 
and the accuracy [64]. To test the ability of the used 
ab initio approach with the employed basis sets and 
energy-consistent pseudopotentials to reproduce the ex¬ 
perimental data, we check the accuracy of the atomic re¬ 
sults. The predicted ionization potentials are 49243 cm“^ 
for Ca, 45814 cm“^ for Sr, 41790 cm“^ for Ba, and 
50073 cm~^ for Yb, in good agreement with the exper¬ 
imental values of 49305 cm~^, 45932 cm“^, 42035 cm~^, 
and 50443 [49] , respectively. The predicted static electric 
dipole polarizabilities of the ground state are 76.5 aj] for 
Ca+, 92.0Og for Sr+, 121.7ag for Ba+, and 61.6Ug for 
Yb+, in good agreement with the experimental values of 
75.3±0.4ag [65], 91.3±0.9a3 [66], 123.88 ± 0.05aj] [67], 



binding energy (MHz) 

FIG. 4. (Color online) Absolute values of the permanent elec¬ 
tric dipole moment of the (GaCr)’'", (SrCr)^, (BaGr)’'", and 
(YbGr)^ molecular ions in different vibrational levels of the 
a®E^ electronic state as a function of the binding energy. The 
origin is in the center of mass. 

and 62.04 flg [47], respectively, and with the most recent 
theoretical results [48]. The potential well depths of the 
alkaline-earth-metal and alkali-metal dimers, using the 
same methodology as in the present study, are repro¬ 
duced with an error of a few percent as compared to 
the experimental results [68, 69]. Based on the above 
and additional convergence analysis, we estimate the to¬ 
tal uncertainty of the calculated potential energy curves 
to be of the order of 100-200 cm“^ that corresponds to 
3-5% of the interaction energy. The uncertainty of other 
electronic properties, including the long-range interaction 
coefficients, is of the same order of magnitude. 


B. Permanent electric dipole moments and static 
electric dipole polarizabilities 

A static electric field, that couples with a permanent 
electric dipole moment and orients molecules, can be used 
to control and manipulate the dynamics of neutral polar 
molecules at ultralow temperatures [70]. Unfortunately, 
the spatial motion of charged ions is highly sensitive to an 
external dc electric field and the control of their dynamics 
with this kind of field is limited. Instead, a nonresonant 
laser field, that couples with a electric dipole polarizbility 
and aligns molecules, can be used to tune ion-atom col¬ 
lision by means of the ac Stark shift. Such control of the 
photodissociation [71], photoassociation [72], magnetoas¬ 
sociation [73], and rovibrational structure dynamics [69] 
has already been investigated in the case of neutral di¬ 
atomics. 

The permanent electric dipole moments of the 
(CaCr)+, (SrCr)+, (BaCr)+, and (YbCr)+ molecular 
ions in the a®E+ and X®S+ electronic states as functions 
of the internuclear distance R are presented in Fig. 3(a) 
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FIG. 5. (Color online) Average static electric dipole polariz¬ 
ability (a) and polarizability anisotropy (b) of the (CaCr)”'", 
(SrCr) + , (BaCr) + , and (YbCr)+ molecular ions in the 
electronic state as a function of the internuclear distance. The 
points indicate the values for the ground rovibrational level. 


and Fig. 3(b), respectively. The values for the ground 
rovibrational level are reported in Table II and the val¬ 
ues for all vibrational levels of the a®E+ electronic state 
are shown in Fig. 4. 

The permanent electric dipole moments are calculated 
with respect to the center of mass, which is a natu¬ 
ral choice for investigating the rovibrational dynamics. 
Their absolute values increase with increasing internu¬ 
clear distance and asymptomatically approach the limit¬ 
ing cases where the charge is completely localized on the 
ion (the dotted lines in Fig. 3). The obtained dipole mo¬ 
ments of the a®S+ and electronic states are sim¬ 

ilar. The difference between the calculated values and 
the limiting cases is the interaction-induced variation of 
the permanent electric dipole moment or, in other words, 
the degree of charge delocalization. This behavior is typ¬ 
ical for heteronuclear molecular ions and implies that 
even molecular ions in very weakly bound states have ef¬ 
fectively a significant permanent electric dipole moment 
in contrast to neutral molecules. In fact, for heteronu¬ 
clear molecular ions the value of the vibrationally aver¬ 
aged permanent electric dipole moment in the body-fixed 
frame with respect to the center of mass increases with 
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FIG. 6. (Color online) Average static electric dipole polariz¬ 
ability (a) and polarizability anisotropy (b) of the (GaCr)"*", 
(SrCr)^, (BaCr)^, and (YbGr)"*" molecular ions in the 
electronic state as a function of the internuclear distance. The 
points indicate the values for the ground rovibrational level. 


decreasing binding energy as shown in Fig. 4 for the in¬ 
vestigated molecular ions in the a®S+ electronic state. 
For levels with binding energies below 1 GHz the electric 
dipole moment exceeds 100 D and reaches 1000 D for the 
most weakly bound levels. Thus, weakly bound molecu¬ 
lar ions have a giant permanent electric dipole moment 
and present excellent prospects for controlling their dy¬ 
namics with the nonresonant laser field. 

The perpendicular a± = a^x = o.yy and parallel 
a II = azz components of the static electric dipole polariz¬ 
ability tensor were obtained in the ab initio calculations. 
The average polarizability a = {2a±_ -|-a||)/3 and the po¬ 
larizability anisotropy Aa = q;|| — a± of the (CaCr)+, 
(SrCr)+, (BaCr)+, and (YbCr)+ molecular ions in the 
a®E+ and electronic states as functions of the in¬ 

ternuclear distance R are presented in Fig. 5 and Fig. 6, 
respectively. The values for the ground rovibrational 
level are reported in Table II. 

The polarizabilities show an overall smooth behavior 
and tend smoothly to their asymptotic atomic values. 
The interaction-induced variations of the polarizability 
are clearly visible while changing the internuclear dis¬ 
tance R and are more pronounced for the elec- 
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tronic state as compared to the a®E+ one. The polar¬ 
izability anisotropies for the rovibrational ground state 
of the investigated molecular ions are slightly smaller 
than the ones of the alkali-metal-atom dimers [74] and 
chromium-closed-shell-atom molecules [36]. Here, we 
present the static polarizabilities which describe the in¬ 
teraction of the molecular ions with the far nonresonant 
laser held. When the shorter-wavelength laser held is ap¬ 
plied the dynamic polarizabilities have to be used, which 
usually are larger but of the same order of magnitude as 
the static ones. 


C. Magnetic Feshbach resonances 

Fermionic alkaline-earth-metal and ytterbium ions, 
and bosonic chromium atoms do not have nuclear spin 
and therefore they do not possess any hyperhne structure, 
in contrast to bosonic alkaline-earth-metal and ytterbium 
ions, and fermionic chromium atom. Three types of mag¬ 
netically tunable Feshbach resonances between the con¬ 
sidered ions and atoms exist, depending on the structure 
of the ion and atom: (i) Feshbach resonances between 
bosonic ions and fermionic atoms are of the same nature 
as those between two alkali-metal atoms [34]. (ii) Fes¬ 
hbach resonances between fermionic (bosonic) ions and 
fermionic (bosonic) atoms result from the interaction of 
the hyperhne structure of the atom (ion) with the elec¬ 
tronic spin of the ion (atom), respectively [9, 59, 60]. 
(iii) Feshbach resonances between fermionic ions and 
bosonic atoms are of the same nature as those between 
two bosonic chromium atoms, that is, they are due to 
the spin-spin interaction [75, 76]. Here, we present ex¬ 
ample calculations for the Feshbach resonances of the 
type (i). The resonances of the types (ii) and (iii) ap¬ 
pear at smaller magnetic held strengths and their typical 
widths are smaller as compared to the type (i) [9, 75]. 

Positions and widths of Feshbach resonances depend 
on the background scattering length, the hyperhne struc¬ 
ture, and the progression of weakly bound rovibrational 
levels just below the atomic threshold. Unfortunately, 
even the most accurate potential energy curves obtained 
in the most sophisticated ab intio calculations do not 
allow one to predict accurately the scattering lengths 
for collisions between many-electron atoms. Thus, at 
present, it is impossible to determine all parameters 
of Feshbach resonances without a priori experimen¬ 
tal knowledge. Nevertheless, the general characteris¬ 
tics of Feshbach resonances such as the density of res¬ 
onances and typical widths can be learned by tuning 
the scattering lengths around the value of the charac¬ 
teristic length scale of the ion-atom interaction. This 
length scale is given by i ?4 = [59, 60] and 

amounts to 1878 bohr for 2212 bohr for 

87gj.-s-^53cr, 2383 bohr for i37Ba+-h53Cr, and 2455 bohr 
for ^^^Yb+d-^^Cr. The characteristic energy scale of 
the ion-atom interaction is given by E 4 = fi^/( 2 /ii? 4 ) 
and amounts to 1.1 /iK for ^^Ca+d-^^Cr, 0.55 /iK for 
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FIG. 7. (Color online) s-wave scattering length as a func¬ 
tion of the magnetic field strength B for collisions of the 
®®Cr (|, —7) atom with the '‘^Ca"'" (4, —3) (a), Sr”*" (5, —4) (b), 
Ba^ (1,1) (c), and Yb'*' (3, —2) (d) ions at the collision energy 
of 100nK. Typical scattering lengths of 0.1i?4, i? 4 , and 2 R 4 
for the electronic state and —0.1J?4, —R 4 , and — 2 R 4 

for the electronic state are assumed, respectively. 


87sr+-p53cr, 0.41/iK for i37Ba+-h53Cr, and 0.37/iK for 
i73Yb+-|-53(^j.^ For collision energies larger than E 4 , 
higher partial waves contribute to the rates for elastic 
scattering and pronounced shape resonances, typical for 
ion-atom collisions, can be expected [9]. 

Figure 7 shows the s-wave scattering lengths as 
functions of the magnetic held strength B for ultra¬ 
cold collisions between the ^^Cr (|,—|) atom and the 
43Ca+ (4,-3), Sr+ (5,-4), Ba+ (1,1), and Yb+ (3,-2) 
ions, respectively. (/, m/) denotes the hyperhne state of 
the atom (ion), where / is the total angular momentum 
of the atom (ion) and m/ is its projection along a quan¬ 
tization axis dehned by the external magnetic held B. 
The results for three sets of the scattering lengths as¬ 
sociated with the and a®S+ electronic states are 

presented. The same absolute values of 0.1i?4, i? 4 , and 




































































































2 i ?4 are assumed for both states with the positive sign 
for and the negative sign for a®!]'*'. The scatter¬ 

ing lengths are set by scaling the potential energy curves 
with appropriate factors A, V{R) XV{R), taking val¬ 
ues in the range of 0.99-1.01. Two families of resonances 
can be distinguished in Fig. 7: broad ones and narrow 
ones. The broad resonances are caused by the hyperfine 
interaction in the ion, whereas the narrow resonances are 
due to the hyperfine interaction in the chromium atom. 
The observed difference in the resonance widths results 
from the difference in the hyperfine coupling constant, 
which is much smaller for the chromium atom as com¬ 
pared to the ions. The widths of the broad resonances 
are of the order of a few to several Gauss. The widths 
of the narrow resonances are relatively small, mostly be¬ 
low 1 G. The density of the narrow resonances is higher 
than the density of the broad ones. Resonances due to 
the dipolar spin-spin interaction, which couples s-wave 
channel to d-wave channel [62], are unnoticeable in the 
scale of the plot. To avoid potential losses due to the 
Zeeman relaxation in an atomic gas [76], the chromium 
atoms in the lowest Zeeman and hyperfine state can be 
used. When choosing the fully spin-stretched hyperfine 
state of the atom (ion), only broad (narrow) Feshbach 
resonances are present in the system. 

Finally, based on the above presented calculations for 
typical scattering lengths and example hyperfine states 
of the ions and atoms, one can expect the existence of 
at least one broad Feshbach resonance and many narrow 
Feshbach resonances at magnetic field strengths below 
1000 G in the case of the investigated type (i) of ion-atom 
mixtures. More Feshbach resonances at smaller magnetic 
field strengths but with smaller widths are expected for 
the types (ii) and (iii) of systems. The probability that 
the true potential energy curves correspond to very unfa¬ 
vorable values of the scattering lengths are rather small 
and such an unfortunate case can be remedied by chang¬ 
ing the hyperfine state or isotope. 


IV. SUMMARY AND CONCLUSIONS 

We have proposed the Ga’*', Sr+, Ba+, and Yb'*' ions 
immersed in an ultracold gas of the Cr atoms as sys¬ 
tems in which ion-atom interactions at ultralow tem¬ 
peratures can be controlled with magnetically tunable 
Feshbach resonances without charge transfer and radia¬ 
tive losses. We have investigated ab initio properties of 
the and a®E+ electronic states of the (CaCr)+, 

(SrGr)+, (BaCr)+, and (YbGr)+ molecular ions. Poten¬ 
tial energy curves, permanent electric dipole moments, 
and static electric dipole polarizabilities have been cal¬ 
culated with the spin restricted open-shell coupled clus¬ 
ter method restricted to single, double, and noniterative 
triple excitations and the multireference configuration in¬ 
teraction method restricted to single and double excita¬ 
tions. Next, magnetically tunable Feshbach resonances 
have been analyzed. 


The ionization potential of the chromium atom is 
larger than the electron affinity of the alkaline-earth- 
metal and ytterbium ions. Therefore, the charge transfer 
process is energetically forbidden and thus the radiative 
losses observed in other experimentally investigated cold 
hybrid ion-atom systems can be avoided in the proposed 
mixtures. The other possible combinations are the Ba"*" 
ions and the Li atoms as well as the Sr+ or Ba+ ions and 
the Dy or Er atoms. The avoidance of radiative charge 
transfer losses opens the way towards robust quantum 
simulations [8] and computations [6] with ultracold het- 
eronuclear ion-atom systems. 

Magnetically tunable Fesbhach resonances have been 
investigated for typical scattering lengths and example 
hyperfine states of the ions and atoms. The assumed 
ultracold collision energies require trapping ions in the 
dipole trap rather than in the Paul trap to avoid heat¬ 
ing due to the rf-field-induced micromotion [1]. The lack 
of the radiative charge transfer losses in the investigated 
systems should enable sympathetic cooling of ions to ul¬ 
tracold temperatures. Presented calculations suggest the 
existence of useful Fesbhach resonances at moderate mag¬ 
netic field strengths but the exact positions and widths 
of these resonances depend on the true values of the scat¬ 
tering lengths which have to be provided by experimental 
measurements. 

The reduction of the two-body losses should also facili¬ 
tate the formation, observation, and application of molec¬ 
ular ions. Magnetically tunable Feshbach resonances can 
be employed in the magnetoassociation [34], which can 
be followed by an optical stabilization to more deeply 
bound rovibrational levels. The photoassociation to ex¬ 
cited electronic states should also be possible [77] but the 
actual pathways of this laser-induced association have to 
be determined. Additionally, the laser excitation of the 
ions or atoms enables radiative charge transfer thus al¬ 
lowing for the study of controlled chemical reactions in 
the investigated systems as well. In the absence of the 
two-body losses, it would be interesting to investigate the 
three-body ion-atom recombination process [78, 79] 

A+-h2Gr^ (XGr)+(u,j)-kGr (5) 

which leads to the formation of the molecular ions 
(ACr)+ in different rovibrational levels (y,j). The con¬ 
trol of this process with magnetically tunable Feshbach 
resonances both between ion and atoms as well as be¬ 
tween atoms should be possible and give new insight into 
few-body physics. One can also think about investigating 
Efimov physics in the ion-atom-atom scenario [80]. 
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